Generation of phase-coherent states 
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An interaction scheme involving nonlinear \ media is suggested for the generation of phase- 
coherent states (PCS). The setup is based on parametric amplification of vacuum followed by up- 
conversion of the resulting twin-beam. The involved nonlinear interactions are studied by the exact 
numerical diagonalization. An experimentally achievable working regime to approximate PCS with 
high conversion rate is given, and the validity of parametric approximation is discussed. 



I. INTRODUCTION 



5-j - 

Optical processes taking place in \ media involve three light waves, yielding to a considerably rich variety 
of nonlinear phenomena, both in the semi-classical jl|] and in the quantum domain j^]. The quantum statistical 
properties of the radiation coming out from such interactions have attracted much attention (see for example Refs. 
||[f |). Squeezing, anti-bunching and entanglement have been predicted and subsequently observed in a series of 
. fascinating experiments [||[|. Most of the theoretical approaches to a quantum theory of three- waves devices have 
' been carried out using the so-called parametric approximation 0. In this framework one of the field modes is in 
\f~} , a strong semi-classical coherent state, so that its depletion as well as its quantum fluctuations can be neglected. A 
fully analytical treatment of the quantum dynamics is not available, whereas numerical methods have been developed 
in the cases of photon number ||] or coherent input states |J. In this paper the three- wave dynamics is evaluated 
without approximation for arbitrary input states, resorting to the numerical block-diagonalization of the Hamiltonian 
in invariant subspaces of the constants of motion. 



In the rotating wave approximation the non-degenerate three-wave interactions are described by the Hamiltonian 

*G ! H oc x (2) [abc) + atht c ] f (!) 

P-i. 

' where a, b and c are the annihilation operators of the three relevant modes, whose frequencies satisfy the relation 
uj c = uj a Depending on the input state of the field, the Hamiltonian (|l]) describes phase-insensitive amplification 

or frequency up- or down-conversion. The first kind of process occurs in situations with small a and b and large 
q-( coherent c, so that the pumping mode can be considered as undepleted and treated as a c-number, in the so-called 
r* ■ parametric approximation. On the other hand, when all the three modes participates to the quantum dynamics we 
. £h ! are m t ne presence of frequency up- and down-conversion processes. 

In the present paper we are interested in the situation depicted in Fig. |l|, where the interaction Hamiltonian (|l|) 
is applied twice: in the first step as a parametric (spontaneous) down-conversion of the vacuum state, generating a 
twin-beam on modes a and b, and in the second step as the up-conversion of twin-beam into mode c. 
This scheme is of interest because, as shown in the following, the outgoing quantum state of radiation closely resembles 
the phase-coherent state (PCS) 

OO 

|A) = v^W E A " I") . ( 2 ) 

n=0 

which has been introduced by Shapiro et al in Ref. Jic[ |. The phase-coherent states are interesting because they are 
optimum phase states for both the Siissmann and the reciprocal peak likelihood |l0[ ] measure of phase uncertainty 
0,^2|. On the other hand, they also could serve as seed state jl^-[l5|| in sampling canonical phase distribution by 
unconventional heterodyne detection p6[ . Moreover, one should mention that the PCS maintain phase coherence 
under phase amplification [ fi"7| ], such that they are privileged states for phase-based communication channels. 
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In suggesting the present scheme we have been inspired by Ref . |Q , where an ideal scheme using a photon number 
duplicator (PND) was suggested for PCS synthesis from twin-beam. As a matter of fact, in such photon recombination 
process the PND is well approximated by the up-conversion from Hamiltonian (EJ) with mode c initially in the vacuum 
jn^Sol. I n this paper we show that the interaction scheme sketched in Fig. EJ is indeed effective for the generation of 
PCS. 

The paper is structured as follows. In Section [n] we briefly describe our approach to the evaluation of the dynamics 
of the three- wave interactions, and discuss the validity of the parametric approximation in the generation of twin-beam 
state. In Section III we analyze the performances of the twin-beam up-conversion in producing phase-coherent states 
in the second stage of the scheme of Fig. El Finally, Section IV closes the paper with some concluding remarks. 



II. DYNAMICS OF THE THREE- WAVE INTERACTIONS 



The overall input state describing the three involved modes can be written in the Fock basis as 

N 

L-Tli ,712 ,713 

\ni,n 2 ,n 3 ) , (3) 

ni,n 2 ,n 3 



where N is an arbitrarily large integer, which denotes the largest not-negligible Fock component. In order to compute 
the dynamical evolution of \?p) 

Wt) = exp (-itH) |^o) , (4) 

one should, in principle, diagonalize the full Hamiltonian matrix in the Fock basis. This becomes a very difficult task 
when the truncation N of the Fock space increases, becoming unrealistic for N exceeding few teens. However, one can 
notice that the Hamiltonian (EJ) admits two independent constants of motion. For the sake of convenience we choose 
them as the following ones 

S = - [a) a + 6 f 6 + 2c f c] , K = a)a + c + c . (5) 

Conservation of S and K means that subspaces corresponding to given eigenvalues of these quantities are invariant 
under the action of the Hamiltonian (EJ), as [H, S] — and [H, K] = 0. In other words, the Hilbert space 'Ha®'Hb®'Hc 
can be decomposed into the direct sum of subspaces that are invariant under the action of the unitary evolution 
operator in Eq. (EJ). Such a decomposition can be written as follows 

H a ®H b ®H c = eJLo ©fc=o ft sfc , (6) 

where 



Hsk = Span { \k — n) ® \s — k — n) ® \n)} 

n G [0, min(k, s — k)] , (7) 

Span{-} denoting the Hilbert space linearly spanned by the orthogonal vectors within the brackets, and \n\) ® \n2) ® 
\n 3 ) = 712,^3} representing the state which is simultaneously eigenvector of the number operator of the three 
modes. The Hamiltonian (EJ) can be consistently rewritten as a follows 

H = Y J h sk , (8) 

sk 

with h s k acting on Ti s k only. Correspondingly, the representation for |^o) can be written in the following fashion 

N s min(k,s—k) 

W=^^ E c k-n,s-k-n,n\k - n,s - k - n,n) , (9) 

s=0 fc=0 n=0 

which enlightens the invariant subspaces structure. In this way one needs to diagonalize the Hamiltonian only inside 
each invariant subspace, thus leading to a considerable saving of resources. 
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As a first application of the above method we look for the conditions under which the parametric approximation is 
justified in describing the process of frequency down-conversion between a strong semi-classical pump and the vacuum. 
The input state is given by |?/>o) = (0,0,0;}, a >> 1 being the amplitude of a considerably excited coherent state. In 
the parametric approximation the pump mode c in the Hamiltonian (Q) is replaced by the c-number a, thus neglecting 
its quantum fluctuations as well its depletion. Within such approximation the dynamics of the input state 1 0, 0) is 
governed by the evolution operator 

L> = exp{Ca t & t -Cab} , (10) 

where £ — —inta, t being the interaction time and k the coupling constant containing the nonlinear susceptibility. 
The evolution governed by Eq. (jlO|) can be easily computed by means of the Baker-Haussdorff-Campbell formula for 



the SU(1, 1) Lie algebra [18 21 2W, and the output is represented by the twin-beam state 

oo 

\X) = V^WJ2 X n \n,n), (11) 



n=0 



where 



X = -itanh(Kt|a|)e #Q . (12) 
In order to check the theoretical results predicted by the parametric approximation we consider the overlap 



o = VW\x) 

between the state g' coming from the exact evolution and the expected twin-beam In Fig. ||(a) we show the 
behavior of the overlap as a function of the scaled interaction time t = Kt for different values of the pump input 
power. In order to evaluate the efficiency of the process we also consider the energy conversion rate ry, which is defined 
as 

= 1 Tr \q' (n a + h b )] 
1 2 Tr [g m h c ] ' V ' 

where g m = \4>o)('4>o\- In Eq. ( |l3|) i] runs between zero and one, the factor 1/2 coming from frequency conversion. 
In Fig. ^(b) we show the behavior of ?y as a function of the scaled interaction time r for different values of the input 
power, as in Fig. ^(a). It is apparent that parametric approximation is valid also for moderate input power, and 
that one has a considerably wide range of values of the interaction time leading to an overlap very close to unit, the 
weaker is the pump, the larger is this range. On the other hand, these values of the interaction time correspond to a 
low conversion rate. By the way, we note that the interaction time leading to maximum conversion rate follows the 
relation T opt oc (n c ) -1 / 3 . 

III. TWIN-BEAM UP CONVERSION 

In this section we analyze the second step of the PCS generation setup reported in Figjl], namely the three-wave 
interaction starting from the twin-beam input state 

oo 

|X> = V^W E X n \n,n,0). (14) 

n=0 



The complex amplitude x 1S confined in the unit circle, and the mean photon number pertaining the state (14) is 
given by 

{XK + n b + h c \ X ) =2| X |7(1-Ix| 2 )- 

The synthesis of the PCS (^) starting from |x) would be easily achieved by having at disposal a device that performs 
the photon number recombination 

\n,n,0) — ► |0,0,n) . 
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Such kind of transformation has been analyzed in Ref. ]!§[] , and has been shown to correspond to the interaction 
Hamiltonian 

H r = (Jb ] {tfb + 1)"5 C + c\tfb + iy?ab . (15) 

Unfortunately, the Hamiltonian (|l5|) cannot be realized by known optical devices. However, one may notice that the 
perfect number recombination 1, 1, 0} —HO, 0, 1) is performed by the Hamiltonian (0), and this suggests to substitute 
the intensity dependent factor in Eq. ( n5[) by i ts expectation value. In spite of this rather crude approximation, the 
trilinear interaction (Q) has been shown |19| , pof to provide a good approximation of the photon recombination in the 
case of a single photon number state at the input. Here, we analyze the case of the input twin-beam state (|l4|). Our 
aim is to demonstrate that the scheme of Fig. [j] is indeed eff ective in sy nthesizing a PCS. As a parameter to evaluate 
the effectiveness of PCS synthesis we use the overlap O = y (X\g out A) between the state 



9 out = Tr afc 



exp 



(-it&) \x)(x\ exp(W)] , (16) 



exiting the \^ crystal in the mode c and a theoretical PCS |A) corresponding to the same mean photon number. In 
order to evaluate the efficiency of the process we also consider the conversion rate 77, defined as follows 

Tr(g°»*n c ) 

V (x\(na+n b )\ X ) ■ 1 } 

In Fig. |3|we show the behavior of the overlap O and the conversion rate rj as a function of the scaled interaction time 
for different intensity of the incoming twin-beam. A remarkable fact is apparent: interaction times corresponding to 
high conversion rate also optimize the overlap between the outgoing state and the theoretical PCS. This means that 
the up-conversion, although only approximated, produces a recombination process which is at the same time efficient 
and quite precise. One should also mention that for the same interaction times one has a small degree of mixing, 
indicating that the outgoing states are quite pure, and minimum reciprocal peak likelihood, thus confirming good 
phase-coherence properties. 

In Fig. |] we report the maximum overlap, along with the corresponding interaction time and conversion rate, as a 
function of the twin-beam input energy N in = (n a +hb). The overlap O slowly decreases versus the input energy N in , 
whereas the conversion rate r\ is almost independent on this quantity, saturating to a value close to 80%. This results 
in a reliable generation of PCS with overlap between 80% and 100%, for outgoing states with energy N ou t — (n c ) up 
to N out — 20 mean photon number. The corresponding reciprocal peak likelihood 6<f> shows the scaling S(f> cx N 01 ^/ 4 
which, though it is worse than the ideal PCS performances S(j) oc N~^ t , it is far superior to the coherent state level 

The interaction time r opt , which corresponds to the maximum overlap, decreases with the input energy N in . By a 
best fit on data in Fig. |we obtained the scaling power-law r op t ~ 1.4A^~ ' 5 . Remarkably, the same scaling is observed 



as a function of the output energy N out , with only a slight change in the proportionality constant r opt ~ 0.9iV ( 



IV. CONCLUSION 



In this paper we have suggested a scheme to generate the phase-coherent states introduced in Ref. |10|. The setup 
involves two nonlinear crystals and it is based on parametric amplification of the vacuum followed by up-conversion 
of the resulting twin-beam, the up-conversion playing the role of an approximate photon number recombination. 

We found that parametric approximation in down-conversion of the vacuum state is valid also for moderate input 
power, and that one has a considerably wide range of values of the interaction time leading to an overlap very close 
to unit, the weaker is the pump, the larger is this range. However, these values of the interaction time correspond 
to a low conversion rate. On the other hand, we found that the up-conversion process involved in the second step of 
the scheme is both power- efficient and quite precise in the generation of PCS. It is a remarkable fact that the range 
of interaction times leading to high conversion rate also optimize the overlap between the outgoing state and the 
theoretical PCS. We have explored the case of twin-beam input photon number ranging from to 54, and we have 
observed a conversion rate about 80%, with an overlap with ideal PCS between 80% and 100%. This corresponds to 
a reliable generation of PCS up to N out — 20 photons at the output. 
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FIG. 1. Scheme of generation of phase- coherent states. 
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FIG. 2. In (a): overlap O = \J (xlff'lx) between the state g coming from the exact evolution and the twin-beam \x) expected 
within the parametric approximation, as a function of the scaled time r = nt for different values of the pump input photon 
number. In (b): energy conversion rate J|asa function of the scaled time for different values of the pump input photon number. 
In both plots different line styles denotes different values of pump intensity: (n c ) = 81 (dot-dot-dashed), (h c ) = 64 (dotted), 
(n c ) = 49 (dot-dashed), (n c ) = 36 (dashed), (n c ) = 16 (solid). The interaction time leading to maximum conversion rate 
follows the relation T op t oc (ri c ) _1 ^ 3 . 
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FIG. 3. Behavior of O and r\ as a function of the scaled interaction time for different intensity of the incoming twin-beam. 
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FIG. 4. The optimized overlap along with the corresponding interaction time and conversion rate as a function of the 
twin-beam input photons. 
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